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Abstract 

In the recently proposed cosmological constant (CC) filter mechanism based on 
modified gravity in the Palatini formalism, gravity in the radiation, matter and late- 
time de Sitter eras is insensitive to energy sources with the equation of state —1. This 
implies that finite vacuum energy shifts from phase transitions are filtered out, too. 
In this work we investigate the CC filter model at very early times. We find that the 
initial big bang singularity is replaced by a cosmic bounce, where the matter energy 
density and the curvature are finite. In a certain case this finiteness can be observed 
already on the algebraic level. 

1 Introduction 

The core of the old CC problem [TJ |2j [3] is the huge hierarchy between the observed 
value Ao ~ 10~ 47 GeV 4 of the CC in the ACDM model and the value A inferred from 
theory. While the latter is difficult to calculate, it is usually expected to be dominated 
by the ultraviolet sector of the theory implying an enormous magnitude |A| 3> Ao- For 
example, the most reliable contribution to A comes from the electro- weak phase transition 
giving rise to the relation |A| ~ (10 2 GeV) 4 ~ 10 55 Ao [3]. A common way to deal with this 
problem is adding a CC counterterm, which subtracts off the huge vacuum energy density 
and permits a reasonable low-energy universe. The price to pay is the enormous finetuning 
of the counterterm, i.e. the precise choice of tens of digits, which is an undesired property 
for any theory. It is important to note that replacing the CC by a dynamical dark energy 
source generally does not avoid the old CC problem [5]. However, it may help modelling 
the late-time accelerated expansion and other cosmological epochs |6|[?ll8l[9"lll0 | lll | ll2 | ll3]. 

Instead of removing the large CC by a static counterterm, we consider a more dynamical 
approach to neutralise the effects of A. It is easy to see that in standard general relativity 
the large CC would dominate the evolution of the universe at early times. Hence, one has 
to introduce nontrivial changes in gravity |14^ [3J [151 E3 or m the matter sector |17| ITHl [19] 
in order to permit a low curvature cosmos. Further ways to address the old CC problem 
can be found in Refs. |20[ l2Tj |22| |23| |24"| |25| 12 6j . where, for instance, quantum gravity 
effects or self-adjustment mechanisms are important. Finally, neutralising the effect of 
vacuum energy was investigated in Refs. |27[ [28. 29, 3Ql|3lJ[32] from different perspectives. 

In this work we continue exploring the CC filter model proposed in Ref. |33j . This 
approach to solve the old CC problem is based on f(R, Q) modified gravity in the Palatini 
formalism |34[ [35j [36l [371 38] , where the gravity action is formulated in terms of the Ricci 



scalar R and the squared Ricci tensor Q [31?| 0U1 EJ S2]- Palatini theories have a more 
involved algebraic structure than the corresponding theories in the metric formalism [43, 
|4"4"1 14"5| l4"6l 14T1 |4"H1 H9J [501 El] > which permits the construction of a filter for vacuum energy. 
Moreover, in our setup we find only second order field equations just as in general relativity, 
thus extra degrees of freedom are avoided. 

So far, only the region where the matter energy density p is smaller in magnitude than 
the CC A has been studied, which includes the radiation, matter and late-time de Sitter 
epochs. We found already on the algebraic level that energy sources with the equation 
of state (EOS) —1 do not contribute to the curvature at leading order, which is just a 
manifestation of the filter effect. In the following we will extend the analysis by discussing 
the very early cosmic epoch where p may be of the order of A or larger. As the main 
result we find that the filter model avoids the occurrence of an initial curvature singularity, 
which is instead replaced by a cosmic bounce with finite matter energy density and finite 
curvature. For a certain case, this finiteness can be seen also on the algebraic level. Finally, 
it is worth mentioning that the absence of an initial singularity has been observed also for 
other Palatini models in Refs. |52[ 153] . and a connection to loop quantum cosmology was 
suggested |54"] . 

The paper is organised as follows: in Sec. [2] we briefly review how to solve the gravity 
field equations in the Palatini formalism, and the procedure will be applied to the CC filter 
model in Sec. |3j We solve the main equations and discuss the cosmological solutions in 
Sec. 0] before we conclude in Sec. 

In this work the speed of light c and the Planck constant h are set to unity, the 
signature of the metric is (—1, +1, +1, +1)- The tensor indices a, ...,/, m, n run from 
to 3, whereas i,j = 1,2,3 are related to spatial coordinates. 



2 f(R, Q) modified gravity in the Palatini formalism 

The CC filter mechanism is based on f(R,Q) modified gravitjQin the Palatini formalism, 
where R is the Ricci scalar and Q the squared Ricci tensor. In terms of the space-time 
metric g a b and its determinant g = det g a b the complete action is given by 



g\\f{R,Q) 



+ S ma t\9abM (!) 



with the matter fields (j) in 5 ma t coupled minimally to g a \,. Furthermore, here the connec- 
tion r^ c and the Ricci tensor R a b are restricted to be symmetric^] 

Rab = ^ab,e ~ ^eb,a + ^ab^fe ~ ^af^eb (^) 

R = g ab Rab, Q = R ab R ab = g ac 9 bd RabRcd. (3) 

In the metric formalism the connection would be identified with the Levi-Civita connection 

Pfec = \g ad {9dc,b + 9bd,c ~ 9bc,d), (4) 

and the Ricci tensor and all derived quantities would be functionals of the metric only. 
Differently, the Palatini formalism treats the connection as an independent object in the 

lr The Palatini gravity model discussed in this work does not introduce extra degrees of freedom. There- 
fore, Ostrogradski-type |55| instabilities do not occur despite the squared Ricci tensor Q in the action. 
2 See Ref. [56j for generalisations. 
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beginning. Therefore, the variational principle 5S = yields two equations of motion for 
gravity, one for the metric, 

fR R m + 2 fQ R m R a ~ 7 = T m \ ( 5 ) 



and another one for the connection, 



g\(f R g mn + 2f Q R mn ) = 0, (6) 



where V a denotes the covariant derivative in terms of the yet unknown connection T^ c 
and fit,Q are t ne partial derivatives of / with respect to R and Q, respectively. 

In the following we briefly summarise how to solve these equations, and we point 
the reader to Refs. |33| 157] for a detailed exposition. Since we focus on the background 
cosmology, the matter sector can be described by a perfect fluid with the energy-stress 
tensor given by 

C = (p+pK/ + (p-A)C (7) 

where p is the energy density and p the pressure of ordinary matter with the 4-velocity 
vector u m normalised by u m u m = — 1. A denotes the large CC vacuum energy density. 

From Eq. ([5j two algebraic equations for the two unknowns R and Q can be derived, 
the first one is the trace of ^ 

f R R + 2Qf Q -2f = T, (8) 

where T = T^ 1 . And the second equation follows from treating Q as a matrix equation 
for the Ricci tensor P := R^ in matrix formj^] 



(2/o) 2 ( P+^I) ■ A'" / + 2/(){p + /)) //„, //" (<)) 
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with 



X 2 := 2f Q (p - A) + f Q f + Y := X - c b ■ yjx* - 2f Q (p + p), (10) 

The quadratic equation ([9]) has the solutions 

c a -2f Q (p+l^i^j =XI + Y^~^, (11) 

where c a ^ = ±1 and the signs of the roots in X and Y follow from requiring that the final 
solution of the model solves Eq. The trace of (11 ) yields 

c a (2f Q R + 2f R ) =4X-Y, (12) 



and at this point we simply eliminate all roots and ±1 factors by squaring Eq. (12) twice, 
((2/ Q fl + 2f R f + 8X 2 + 2f Q (p + p)f = 36X 2 (2f Q R + 2f R f . (13) 



Using ([8|) and (13), the Ricci scalar R and the squared Ricci tensor Q can be determined 
as algebraic functions of p, p and A. 



3 Here, / = is the identity matrix, and P 2 = R^R". 
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Next, we derive the connection T with the help of the auxiliary metric h mn defined by 

VW\h mn = V¥\g ma (fJ + VQp) n a , (i4) 

which allows writing Eq. ^ in the form 



V a [T] y/\h\h 



0. 



(15) 



Since P follows from Eq. (11), one finds after some algebra that h mn is related to the 



"physical" metric g mn by the disformal transformation [57 

L 2 



h. 



n I g 



Li - L 



UmUn 



Li 



where 



fi:= VlLiCLi-La)!, 



Li := c a X + -/r, L 2 := c a Y. 



(16) 



(17) 



Finally, the Palatini connection V is just the solution of ( 15 ), i.e. the Levi-Civita connection 

of hm/n. 



l 



'FbcM — ~h a (hdc,b + hbd,c — hb Ci d) 



(18) 



with hmn given in (16). This result solves Eq. ([6| and eventually defines the Ricci tensor 
and all quantities derived from it. Note that h mn does not contain space-time derivatives 
because it is a function of only g a b, p, p and A. Hence, R a bi R an d the other quantities 
calculated according to Eqs. (J2j and ([3| contain derivatives of at most second order of the 
metric g a b- Consequently, in this modified gravity theory in the Palatini formalism one has 
to solve the second order differential equation 



R[T]=R( P , P ,A), 



(19) 



where the right-hand side follows from the algebraic solution of Eqs. (|8|) and (13), whereas 



the left-hand side is given in Eq. ([3]) using T from Eq. (18). In the next section, we will 
apply this procedure to the CC filter model. 



3 CC filter model 

In Ref. [33j a modified gravity model was constructed that filters out vacuum energy 
contributions (i.e. everything with the EOS p = —p) to the CC. It is defined by the 
following action functional 

( 2 \ m 
-^J , B:=R 2 -Q, (20) 

with the dimensionful parameters k and (3 and the positive number m. Note that in this 
setup, the term z is never considered to be a little correction to kR, instead both terms 
in / act in collaboration. From here one, we will use B instead of Q. With the partial 
derivatives of /, fu = (kR+^itiz) / R — 2/qR, /q = mz/B, and the trace T = — 4A + 3p — p 
of the energy-stress tensor ([7| , Eq. ^ becomes 

4 

7z = kR- 4A + 3p- p, j:=-2--m, (21) 
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while Eq. ( 13 ) can be written as 







2 

= kR + r mz 

9 



+ 



27 mZ 



'FT 



\R 2 ) 
i + 



i + 

3kR 



— (3kR + 2r) + T7 9 , 9 ((kR) 2 - r 2 
mz 4(mz) 2 y 

2 



(22) 



2{mz) 



In the last equation we introduced the variable r := p + p, which vanishes identically for 
CC contributions. 

For matter energy densities p and curvatures kR smaller in magnitude than the large 
CC, the term "fz ~ — 4A in (21) is approximately constant and Eq. (22) determines the 



curvature as a function of r alone. Therefore, contributions to the CC and shifts in the 
vacuum energy density are filtered out. This limit of constant z = 0(A) covers most of the 
radiation epoch, the matter era and the late-time de Sitter phase, and it has been analysed 
analytically in Ref. |33| by neglecting strongly suppressed terms like 0(r/z). 

In the following we complement our previous results by investigating the very early 



universe, where r may be of the same order as A or above. Consequently, z in Eq. (21) is 



not constant anymore and all terms in Eq. (22 ) must be taken into account. It is convenient 



to work with dimensionless variables, where energy densities and the curvature term kR 
are normalised by the cosmological constant (4A), 



kR r 
y:= 4A' X:= 4A 



Moreover, let us define 



m ■ z 



m 



m. 



[l-y-x-u) n 



P + P 
4A ' 



with u r , 



3co - 1 



4A —7 + 1 

where lv = p/p denotes the matter EOS. Finally, we introduce the variable 

R 2 



W :-- 



B 



ya (1 - y - x ■ u m ) m • 5, 



(23) 



(24) 



(25) 



where B was eliminated on the right-hand side by using the definition of z in Eq. (20) 



The constant 5 is completely fixed by the model parameters and the large vacuum energy 
density A, 



I 2l " I f ' ) , where (p e )s 



9 77 V Pe J 



4 (2 4A\ /-7/3 
k3 -m — — — — 
\9 -jj V 4A 



(26) 



According to our earlier results in Ref. [33J, 5 is very large in magnitude because p e is 
of the order of the late-time critical energy density in the limit 72; ~ — 4A. Expressed in 



terms of the dimensionless variables the main equation (|22|) reads 

8 m. 







y + x 



2 m z 
9VV 



i + 3(3y + 2s) 9(y 2 



x 



2M 



m z Ami 
For the late-time solutions mentioned above we find m z 



+ 



27 W 2 



1 + 



3y_ 

2m 7 



-1 2 



(27) 



\W\ 3> 1. Therefore, Eq. (27) simplifies considerably. 

y + x 



2 m z 
9 W 



0. 



-m/7 = O(l), \y\, \x\ 1 and 

(28) 



5 



indicating kR — > p e for vanishing matter r — > 0, which justifies introducing p e in Eq. (26). 



In the next section, we solve the complete equation (27) requiring a numerical treat- 
ment. In practise, we work along the following procedure. First, we consider a spatially flat 
cosmo logical background with the expansion described by the scale factor a(t) as a func- 
tion of cosmological time. In Cartesian coordinates (i = 1,2,3) the non-zero components 
of the metric are given by goo = — 1> 9u = a2 (t)i an d m the stress tensor we use u m = 5^ 



. Hence, the auxiliary metric in ( 16 ) has only the diagonal components 



oo 



n 



L 



12 



hi 



fla 2 (t), Lu := Li - L 2 , 



(29) 



which determine via ( 18 ) the symmetric components of the Palatini connection 

L[(x) L' 12 (xY 



1 oo 



r 



X 

4 

a* 
Li 
a 



Li 

a T x 
a 4 



;o 



X 

a + 4 



12 



, X 



L\2 

L' 12 (x) + 
L[(x) 



L 12 L[(x) 
Li 



L 



+ 



12 



(30) 



Here, the overdot indicates a derivative with respect to the time coordinate t, whereas the 
prime denotes the partial derivative with respect to x. For simplicity we consider matter 
with a constant EOS u, which obeys the standard conservation equation p + 3^p(l+u) = 
implying p oc a -3 ( 1+aJ ) and x = — 3-(l + uj)x via x = (1 + cj)/9/(4A). Consequently, the 
Ricci scalar R(t) calculated via (J3| from the connection (30) is a function of a, a and a, 
whereas the algebraic solution y{x) = kR/{AK) of Eq. (27) is only a function of a. For 



obtaining a(t) we solve numerically the differential equation (19) 



m 

H 2 



K-y(x(t)) with K 



4A 



(31) 

wc 
In 



where Hq = ^(^o) is the Hubble rate at the arbitrary initial time to- in Ref< |3 
found (— k) to be of the order of the inverse Newton's constant G indicating kHq 
by identifying Hq with the present Hubble rate and using the late-time solution ( p8_ 
this case, K is roughly the square root of 5 in (26), its exact value is fixed when we 
specify Hq and x(io) at the initiaQ time to. Numerically, we first solve the equation^] 
R 2 /B = W(y(x)) for a(to)j where R 2 /B has to be calculated from (30) at t = to an d 
x = x(ao). Then, with a(to) the values of R(to) and y(x(to)) in Eq. (31 ) are known, which 



yields the numerical value of K and Eq. (31) can be solved to determine the evolution of 
the scale factor a(t) for all times t. 



4 Solutions 



Here, we study the cosmological evolution resulting from the differential equation (31). 



Thereby we use the solutions to the main equation (27), which relate the curvature y 



Alternatively, one could start with given values of k and x(to) and infer Ho from it. This is analogous 
to finding Ho from the Friedmann equation 3Ho = 8irG ■ p{to) in general relativity. 



3 We use this equation in terms of arbitrary t to determine the correct choice of c a ,b in (10 1 and (171 



after Eq. (31 1 has been solved for a(t). 



I. 




to the matter energy density and pressure in x as given in (23). We do not restrict the 



magnitudes of x and y, however, in order to be roughly compatible with our existent 
universe, we consider only solutions y(x) with both a final accelerating phase, x — > 0, and 
a high energy decelerating epoch, y ~ —x. Moreover, for the numerical examples below we 
set m = 1/3 in Eq. (20), because this parameter has only little influence on the qualitative 
behaviour of the solutions. With these preliminaries, there are four qualitatively different 
cases to study corresponding to A gs and 5 gs 0. Note that the sign of 5 can be adjusted 



by a suitable choice of the parameter /3 in the action functional ( 20 ) . 



1. The absence of an initial singularity can be seen best in the case 5 > and A < 



implying p e < via Eq. (26) and 7 < from (21 ). In this case the matter variable x 



r/(4A) is negative because of r = p + p > 0. Numerically, we do not find solutions 
for y < 0, that satisfy our requirements given above. This can be inferred also from 



Eq. (28 ), where all terms would be negative. On the other hand, for y > (or kR < 0) 
we find Eq. (28) as a good low-energy description with y(x — > 0) — > p e /(4A) and 
y ~ —x for \p e \ <C \x\ <C 1. However, before y — > 1, the matter variable |x| becomes 
smaller again and eventually vanishes for finite y < 1 as illustrated in Fig. [I] As a 
result, \x\ oc a _3 ^ +1 ^ and y are bounded from above signalling the absence of the 
initial singularity, where x or y would diverge. Analytically, this boundedness can 



be understood from Eqs. (24) and (25) indicating that for sufficiently large values 
of y or |x| the variables W and m z vanish. Since they appear in the denominators in 
our main equation (27), it is clear that neither \x\ nor y can become larger than 1. 



After solving the differential equation ( |31[ ) for the scale factor a(t), we observe that 
not the whole phase-space curve y(x) in Fig. [I] is realised, but only a smaller part of 
it. This behaviour can be explained by observing that the function L\ — L2 vanishes 
at x ~ —0, 043 as shown in Fig. [2j Since L\ — L2 appears in the denominator of the 



connection components (30) and therefore in the Ricci scalar R(t), the equation (31 ) 
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Figure 2: The functions L\{x) and L\2{x) in the connection components in Eqs. (30) for 
the case 5 = 5- 10 6 , A < and two values of the EOS u> = g, 0. Notice the zero of L\2 at 
x ~ — 0, 043, where the cosmic bounce occurs. 



can be fulfilled only if a/a oc x/x vanishes at the zero of L\ — L2 because y(x(t)) 



on the right-hand side of (31) is always finite. Accordingly a = indicates the 
occurrence of a cosmic bounce, which is shown in the plots for the scale factor a(t), 
the Hubble rate H = a/a and the deceleration q(t) in Fig. [3} Also, we can confirm the 
analytical results from Ref. |33j . where it was found that at early times the universe 
behaves like a radiation- or dust-dominated cosmos depending on the matter EOS w. 
Thus for uj = g or the deceleration q approaches the value 1 or |, respectively. In 
this era, y ~ —x is a good approximation. Finally at late times (x — > 0) we find a 
de Sitter cosmos with y — > p e /(4A). 

2. The second case we consider is 5 < and A > implying x > and p e < 0. Here, 
the only reasonable solution requires y, kR < 0. In contrast to the first case, the 



variables x and \y\ are not bounded from above because W and m z in Eqs. (24) 



and (25) do not vanish. Instead we find y — > —x for large x — > 00, see Fig. [I] Also, 



the functions L\ and L\ — L2 do not vanish in this region. Nevertheless, the cosmic 



evolution inferred from the numerical solution of (31) shows a cosmic bounce, where 
the scale factor a(t) and the matter and curvature variables x,\y\ remain finite, see 
Fig. |5} As before, the universe shows a radiation- or matter-dominated expansion 
behaviour with y ~ — x and a final de Sitter phase (y — > p e /{AK)). 

3. The third possibility is characterised by 5 < and A < indicating x < and p e > 0. 
We did not find a reasonable solution in this case, although for y > a y ~ —x phase 
exists, but it cannot end in a low-energy accelerating phase with x — > because all 



terms in Eq. (27) are positive in this limit. Hence, we do not pursue this solution 
any further. 

4. Finally the last case with 5 > and A > implies x > 0, and for y < we 
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Figure 3: The numerical solutions from Eq. (31) for the scale factor a(t), the Hubble 
rate H(t) and the deceleration factor q = — |/IPfor the case 5 = 5- 10 6 , A < described 
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Figure 4: The solution y(x) from Eq. (27) and the functions L\(x) and L\2{x) in the 
connection in Eqs. (30 ) for the case 8 = — 5-10 6 , A > and two values of the EOS u = \,Q. 
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find a radiation/matter phase described by y ~ —x for x — > oo. Similar to the 
third case, p e > does not seem to allow a late-time de Sitter solution, however, 



the last term oc W~ in Eq. (27) is positive giving rise to an accelerating solution 
with x — > 0. This cosmos is, however, quite exotic, because there is a contraction 
phase between the decelerating and the final de Sitter phase. For completeness the 
expansion behaviour of this interesting but unrealistic solution is shown in Fig. [6] 

In the examples shown in this section we used relatively small values of |<5| for numerical 



reasons. According to its definition in Eq. (26), 5 ~ (A/p e ) 7 ^ 3 is related to the ratio of the 



large CC and the late-time energy density implying an enormous magnitude for 6. However, 



this does not correspond to unreasonable values for the parameter f3 in the action (20), 
instead a range of common values can be obtained depending on m and k. Moreover, in 
this work we have normalised energy densities like r and kR by the constant A such that 
our results are independent of the numerical values of k and A. For details on (5 and the 
relation between k and Newton's constant we refer the reader to Ref. [33J. 

In summary, we found cosmological solutions for both signs of the large CC A, where 
the initial big bang singularity is replaced by a cosmic bounce with finite curvature and 
finite matter energy density. After the bounce these universes behave like being radiation 
or matter dominated, where the deceleration q is close to 1 or ^, respectively. Eventually, 
when the matter component is sufficiently diluted the final de Sitter era begins, with an 
effective vacuum energy density ~ p e much smaller in magnitude than the large CC A. 
The latter property and the absence of an initial cosmological singularity are the main 
benefits of the CC filter model. 
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Figure 6: The solution x(y) from Eq. (27), the scale factor a(t) and the Hubble rate H(t) 



from the numerical solution of Eq. (31) for the exotic case 5 = 5000, A > 0. Notice the 



contraction phase between the decelerating and accelerating eras. 



5 Conclusions 

The filter model for a large CC in the context of Palatini modified gravity |33| makes 
gravity insensitive to vacuum energy density contributions. In other words, the curvature 
is dominated by ordinary matter, whereas energy sources with the EOS of a CC do not 
contribute at leading order despite having a much larger energy density at late times. 

In this work, we have explored the CC filter model in the context of the very early uni- 
verse, where in classical theories one usually expects to find an initial big bang singularity 
at very high energies. However, in our setup we observed that the curvature singularity 
is replaced by a cosmic bounce with the matter energy density, the scale factor and the 
curvature remaining finite. In the case shown in Fig. [T]the boundedness can be seen already 
from the solution of the algebraic equation ( |27[ ), whereas the other cases show the same 
feature after solving the differential equation (|31[). Moreover, our numerical results also 
confirm the analytical study of the late-time evolution in Ref. [33J. 

The absence of an initial space-time singularity could mean that our classical model 
for gravity does not require an ultraviolet completion from a theory of quantum gravity. 



However, in the light of the exotic form of the gravity action functional (20) it is unlikely 
that this setup is a fundamental theory. On the other hand, we may follow the argu- 
mentation in Ref. |54j . where Palatini gravity models are proposed as effective continuum 
descriptions of loop quantum cosmology. In both approaches the regularisation of the ini- 
tial singularity does not require the introduction of new degrees of freedom making them 
minimal solutions to this problem. 

Independent of its theoretical origin, the CC filter model has the capability to handle 
two fundamental problems of gravity, the old CC problem and the avoidance of an initial 
singularity. Therefore, it is an attractive scenario for further investigations. In addition 
to the qualitative features discussed so far, it will be interesting to study more aspects 
of the CC filter model in the future, e.g. black hole or star interiors or anisotropics and 
inhomogeneities in cosmology. 
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